Numerical Results for U(l) Gauge Theory on 
2d and 4d Non-Commutative Spaces 



W. Bietenholz l , A. Bigarini 2 , F. Hofheinz 
J. Nishimura 3 , Y. Susaki 3 4 and J. Volkholz 1 * 

1 Insistut fur Physik, Humboldt Universitat zu Berlin 
Newtonstr. 15, D-12489 Berlin, Germany 
2 Dipartimento di Fisica, Universita' degli Studi di Perugia, 
and INFN, Sezione di Perugia, Via Pascoli 1, 1-06100 Italy 
3 High Energy Accelerator Research Organization (KEK) 
1-1 Oho, Tsukuba 305-0801, Japan 
4 Institute of Physics, University of Tsukuba, 
Tsukuba, Ibaraki 305-8571, Japan 



Abstract: We present non-perturbative results for U(l) gauge theory in spaces, which 
include a non-commutative plane. In contrast to the commutative space, such gauge the- 
ories involve a Yang-Mills term, and the Wilson loop is complex on the non-perturbative 
level. We first consider the 2d case: small Wilson loops follows an area law, whereas for 
large Wilson loops the complex phase rises linearly with the area. In four dimensions 
the behavior is qualitatively similar for loops in the non-commutative plane, whereas the 
loops in other planes follow closely the commutative pattern. In d = 2 our results can 
be extrapolated safely to the continuum limit, and in d = 4 we report on recent progress 
towards this goal. 



1 Non- commutative U(l) gauge theory 

We consider the simplest version of non-commutative (NC) spaces: only two (Euclidean) 
coordinates are NC, and the non-commutativity parameter 9 is constant, 

=i0€ij (indices i,j G {1,2}) . (1) 

According to a historic observation by Peierls P, such coordinates describe a charged 
particle moving in a (commutative) plane, which is crossed by a strong, orthogonal mag- 
netic field B. In fact, if we ignore the kinetic term in the Lagrangian, such a particle has 
the canonical momentum pi = qBeijXj, where q is the electric charge. Under canonical 
quantization we obtain 

[xi,pi] =ih = qBeij\xi,x 3 \ , (2) 

which is a NC relation with 6 oc 1/B. Based on this idea, NC field theory plays a role as 
a formalism in solid state physics. 
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A similar idea is also used to map string theory in a magnetic background onto NC 
field theory. The construction of this mapping 2. has tremendously boosted the interest 
in NC field theories (i.e. field theories on NC spaces). 

An important mathematical result states that we can return to the use of ordinary 
(commutative) coordinates if all the fields are multiplied by star products, 

<j)(x) -k ip(x) := <j)(x) exp ( - *d i 9 e„ dj ) i/j(x) . (3) 

Here we focus on pure U(l) gauge theory, which has the Euclidean action 

S[A] = \J tfxF^KF^ , 

F» v = d^A v - 8 V A^ + ig[A m A v ]t , (4) 

where the last term is a star-commutator. This action is star-gauge invariant, i.e. invariant 
under transformations 

A„(x) -> U(x) * A^x) * U( X y - -U(x) * dJJ{xy , (5) 

where U(x) is star-unitary, U{xy * U(x) = 1. 

Other U(N) gauge theories may be studied along the same lines, but the formulation 
of SU(N) gauge theories runs into trouble on NC spaces. To see this problem, let us 
represent the gauge field as A^T a , T a being the Hermitian generators of the gauge group. 
Then the star commutator can be decomposed as 

2 [A°T a , A%1 = [Al, AH ■ {T a , T b } + {A% A% • [T a , T h ] . (6) 

If we deal with SU(N) gauge theory we have the additional condition TrT a = 0, which 
should be reproduced by the star-commutator for the algebra to close. In the commutative 
case, 9 = 0, the first term vanishes ( [A a , A^] = ). Then the algebra closes because 
Tr [T a ,Tb] trivially vanishes. For finite 9, however, also the first term contributes, which 
causes trouble (because Tr {T a ,Tb} oc 5 a b does not vanish). 

Therefore it is motivated to concentrate on U(l) as a physical gauge group, which can 
be accommodated on a NC space. 

2 Mapping onto a twisted Eguchi-Kawai model 

As a first step towards a formulation to be used in Monte Carlo simulations, we introduce 
a lattice structure in the NC plane. Since all space-points are somewhat fuzzy, we cannot 
expect sharp lattice sites, but the operator identity 

exp (^i — x^j = 1 (7) 

does nevertheless impose a lattice structure with lattice spacing a. If we require the 
momentum components hi to be periodic over the Brillouin zone, the above condition 
implies that only discrete momenta occur (see e.g. 0), which is characteristic for a finite 
volume. If we assume a periodic N x N lattice, the allowed momenta are separated by 
Aki = -tjT, and the non-commutativity parameter can be identified as 



9 = -No 2 



(8) 



Therefore we are interested in a double scaling limit a — > 0, iV — > oo with iVa 2 = const., 
which leads to a continuous non-commutative plane of infinite extent. 

However, even on the lattice it is far from obvious how to simulate NC gauge theory; 
note that this seems to require star-unitary link variables. It is highly profitable to map 
the system (or its NC part) onto a twisted Eguchi-Kawai model. The latter is defined on 
a single space-point and its action has the form jl] 



U\ and U2 are unitary N x N matrices which encode the degrees of freedom of the lattice 
model, and the twist factor is chosen here as Z21 = Z* 2 = exp(iri(N + 1)/N), where N 
has to be odd. Then there is an exact equivalence to the lattice NC gauge theory, i.e. the 
algebras are fully identical, as Ref. jHj showed in the large N limit. A refined consideration 
found such a mapping even at finite N [6j. This is the form which is suitable for numerical 
simulations. 

It is straightforward to formulate (the analogue of) Wilson loops in the framework of 
the matrix model, 



This corresponds to a rectangular Wilson loop of sides al and aJ . Mapping this quantity 
back to the lattice leads in fact to a sensible definition of a Wilson loop in the NC gauge 
theory jjj. The exponent of the twist factor is taken positive (negative) for an anti- 
clockwise (clockwise) orientation. The Wilson loop is complex in general, although this 
property cannot be seen perturbatively in d = 2 jHJ. For further perturbative studies we 
refer to Refs. 0, and to Refs. jTU| for semi-classical approaches. 

3 Numerical results 
3.1 Two dimensions 

In the planar limit of the two dimensional model we obtain the commutative U(N — > 00) 
lattice gauge theory on a plane, which was solved by Gross and Witten by studying 
Wilson loops [TT]. In this limit they found an exact area law. We tested that the values 
of N, which are accessible to our simulations, do approximate the planar limit well by 
checking the validity of the large TV Schwinger-Dyson equations, along with the large N 
factorization. On the lattice, they relate Wilson loops of different shapes in the planar 
limit. An example for the corresponding contours is illustrated in Figure ^ This was the 
property that the original work by Eguchi and Kawai (without twist) was based on [T2] . 
Indeed, we observe that our measurements for the two sides of such equations converge 
as we increase N to a magnitude of a few hundred, keeping (3 fixed, for instance at 0.25. 
This is shown in Figure 121 see also Ref. [T3] . 

Hence we can use this limit to identify the lattice spacing as (3 oc 1 /a 2 , as in the result 
by Gross and Witten. Therefore we take the double scaling limit by keeping the ratio 
N/P constant. 

In Figure H3 we show our simulation results ^3] for the Wilson loop 




(9) 




(10) 



W(I) = (W 12 (I x /)) 



(11) 




Figure 1: An example for a set of contours, which are involved in a Schwinger-Dyson equation. These 
equations relate the corresponding Wilson loops in the planar limit. 
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Figure 2: The convergence towards the validity of the Schwinger-Dyson equations as N increases at fixed 
(3 = 0.25. On the x-axis of the plot on the left we show the area of the contour called C in Figured On 
the right we show the convergence in N at a fixed area of 7 x 8 for the contour C . Note that in a common 
normalization — with a pref actor 1/N for the Wilson loops — our result is consistent with the expected 
finite N artifacts of the order of 1 / N 2 . 

in polar coordinates, as a function of the physical (i.e. dimensionful) area (al) 2 . We 
choose N over a broad range, with a fixed ratio N/ (3 = 32. 
From Figure El we infer the following properties: 

• The observable W(I) does indeed stabilize in the double scaling limit. Thus we 
have found a new universality class. In particular this shows that the model is 
non-perturbatively renormalizable. 

• At small area, the absolute value follows the Gross- Witten area law, which is marked 
by a line in the plot on the left-hand-side of Figure El In that regime the phase is 
practically zero, hence Wilson loops with areas A = (al) 2 <4 follow the behavior 
the behavior of commutative large iV gauge theory. 

• For larger areas, the absolute value does not decay any further, but the phase starts 
to increase linearly in the area. 

Regarding the second point — the area law for small Wilson loops — we also measured 
the Creutz ratio, 

X(J ' J) = - l0g W-l)xJ)^(/x(J-l)) ) ' (12) 

which singles out the string tension a for decays oc exp(— a A). Typical results for (nearly) 
square shaped Wilson loops, x{I, I)> as we U as extremely anisotropic (rectangular) Wilson 
loops, x(2, J), are shown in Figure |U on the left. We see a stable value a ~ 1 as long 




Figure 3: We show square shaped Wilson loops W(I) on a NC plane: from the plot on the left we see 
that small loops follow the Gross- Witten area law (which is marked by a line). At large areas the absolute 
value does not decay any further, but a linearly rising phase sets in, as the plot on the right shows. 



as we are well inside the area law regime, and a marked deviation from it as the area 
approaches its critical value. 

In view of the third observation — the linearly rising phase — we performed further 
tests with different ratios N/f3 = 16, 24, 48, i.e. different values of 9, see Figure |U on 
the right, and also with rectangular loops. At large areas we always found to a very high 
precision the simple relation 

Phase — — — A - B . (13) 


In the last term we introduced a formal magnetic field B = 1/9, in the spirit of Peierls' 
map. 
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Figure 4: On the left: the Creutz ratio for Wilson loops with small areas, up to the transition observed in 
Figure^ On the right: the linearly increasing phase at various values of the non-commutativity parameter 
9, in agreement with eq. \T$) . 



This behavior corresponds exactly to the Aharonov-B ohm effect. It is known from 
the corresponding magnetic fields in string theory and in solid state physics, where it is 
implemented on tree level. Here, however, we recover the very same behavior unexpectedly 
clS cL dynamical effect at low energy. 

One may wonder why the short-ranged non-commutativity has striking effects on the 
large rather than the small Wilson loops. Apparently there is some UV/IR mixing J3| 
going on, even though there are no divergences in the perturbative expansion of this model. 
This suggests that UV/IR mixing occurs non-perturbatively, and it belongs therefore to 
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Figure 5: On the left: The real part of the plaquettes, averaged over all possible orientations. The dashed 
lines show the asymptotic behavior expected from expansion in strong resp. weak coupling. On the right: 
The real part of the plaquettes, plotted separately according to their orientation. 



the fundamental nature of NC field theory. This is in agreement with analytic ^Bj and 
numerical results for the NC A0 4 model. 

3.2 Four dimensions 

We now proceed to four dimensions, which include a NC plane, as well as two commutative 
coordinates — among them is the Euclidean time. Therefore we deal with one NC plane, 
one commutative plane and 4 mixed planes, and we will distinguish the (flat) Wilson 
loops accordingly. 

There are 1-loop studies about the photon dispersion relation in such spaces [TH] . 
Computations of the (3 function indicate asymptotic freedom [TH]. Other perturbative 
studies suggest that this model may have a negative IR divergence, which would render 
it unstable and unsuitable to describe physics (2Hj. More precisely, those works suggest 
that bosonic (fermionic) degrees of freedom contribute negative (positive) IR poles, and 
they therefore proceed to SUSY models. Here we stay in the simple bosonic framework; 
of course we would like to verify this claim on a non-perturbative level. 

Again we introduce a lattice formulation and map the NC plane on a twisted Eguchi- 
Kawai model; one on each lattice point of the remaining (commutative) plane. On the 
technical side, it is very important to apply the heat bath algorithm for the simulations. 
This is not straightforward, since this matrix model is non-linear in the link variables. 
For our 2d model a method to linearize it with the help of auxiliary fields was introduced 
by Fabricius and Haan 21j. We also applied it in our simulations, which led to the results 
in the previous Subsection. In d = 4 it turned out that this method is still applicable in a 
generalized form [22], which allows us to use the powerful heat bath algorithm. Without 
it, it is very hard to overcome the problems related to thermalization and autocorrelation 

We first present simulation results at iV = 25 j22j. 1 The left side of Figure El shows the 
real part of the plaquette values at different (3. In this case all six possible orientations of 
the plaquettes are averaged over. The results match the anticipated asymptotic behavior 
jl] of the plaquette expectation value (P): at weak coupling is behaves as ~ 1 — and 
at strong coupling as ~ (3. A phase transition occurs around (3 ~ 0.35. All the results 
that we show in the following were obtained in the weak coupling phase. 

1 More precisely, this means that the lattice size is 25 X 25 in the NC plane (remember that this extent has to be odd) 
and 24 X 24 in the commutative plane (this is suitable for parallel computing). 
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Figure 6: Temporal correlation of the plaquettes lying completely in the NC plane (on the left) and 
completely in the commutative plane (on the right). In the former case only we see a trend towards an 
exponential decay. 



The decay of square shaped Wilson loops with increasing area is shown in Figure El on 
the right. In that plot the results are separated according to the plaquette orientations 
(NC, mixed and commutative). The real part of the NC Wilson loops performs a damped 
oscillation around zero. This is in qualitative agreement with the observations in d — 2. 
For the other two types of Wilson loops the real part decays monotonously towards zero. 

We have again measured the Creutz ratios for various rectangular Wilson loops [24 . 
The (dimensionless) string tension approaches zero as (3 —>■ oo. 

Figure El (on the left) shows the correlation function of the plaquettes lying completely 
in the NC plane, separated by At in Euclidean time. These data were taken at /3 = 0.4, 
i.e. barely in the weak coupling phase. The decay is likely to be exponential, but a higher 
statistics is required to verify this behavior. 

The analogous plot for plaquettes in the commutative plane is shown in Figure HI on 
the right. In that case the decay does not seem exponential. 

At last, we illustrate our progress in the search for a double scaling limit. In contrast to 
d — 2, we do not have an analytic result in the planar limit to identify the physical lattice 
spacing (i.e. to relate (3 and a). However, we observed that double scaling is approximated 
quite well if we assume ad hoc the simple relation (3 oc 1/a, so that the double scaling 
limit keeps N/(3 2 = const. In particular we set this ratio to 20 and obtained the results 
for the real parts of square shaped Wilson loops shown in Figure [71 We are also working 
on a systematic search for double scaling by matching the data without any assumption 
about the relation between (3 and a. This represents a completely unbiased test of the 
result presented here. 



4 Conclusions 

We have simulated NC U(l) field theory in d = 2 and in d — 4. In d = 2 we observed an 
area law for Wilson loops with small expected from the equivalence to commu- 

tative large N gauge theory in the UV regime. At large phase starts to increase 

linearly in the area, which is reminiscent of the Aharonov Bohm effect. 

The simulations in d = 4 (with two NC coordinates) are still on-going. We first 
measured the action, plaquettes, Wilson loops and Creutz ratios. The results for the 
action and plaquettes agree with asymptotic predictions. There seems to be a vanishing 
string tension at weak coupling in the part of the parameter space we explored. The 
Wilson loops in the NC plane show a behavior similar to the 2d NC model, whereas the 
loops in other planes are close to the commutative feature. 
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Figure 7: We show the real part of Wilson loops in the three types of planes ( commutative, mixed and 
NC) of our 4d space. Here we have fixed N//3 2 = 20. The results approximate a double scaling behavior 
decently well. (We always refer to N x N matrices for the NC plane, and a (N — 1) x (N — 1) lattice on 
the commutative plane.) 

We found a simple ad hoc ansatz which seems to provide quite well a double scaling 
limit. This ansatz is still being checked with a fully systematic approach. The double 
scaling limit corresponds to the model at zero lattice spacing and in infinite volume, hence 
its identification is a milestone in the non-perturbative analysis of this model. Thus we 
hope to glimpse at the dispersion relations for the NC photon, which might allow us to 
confront the theory with phenomenological data. 

Even before that, however, it is important to clarify if this model is actually IR stable. 
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